Z9  t  l  8  t  Vcw 


At>F366"*0O 

|ad  ~ 

TECHNICAL  REPORT  ARBRL-TR-02496 

A  SURVEY  OF  PENETRATION  MECHANICS 
FOR  LONG  RODS 

Thomas  W.  Wright 


ILi  j 

C 


88  07  26  065 


Destroy  this  report  when  it  is  no  longer  needed. 
Do  not  return  it  to  the  originator. 


Additional  copies  of  this  report  may  be  obtained 
from  the  National  Technical  Information  Service, 

U.  S.  Department  of  Commerce,  Springfield,  Virginia 

22161. 


The  findings  in  this  report  are  not  to  be  construed  as 
an  official  Department  of  the  Army  position  unless 
so  designated  by  other  authorized  documents. 


TJm  uo*  trad i  now*  or  •namf a« r«  • 
*•*  aoniHtkU  i<\aor*mmnt  of  any  a 


name  in  thio  report 
rarial  produot. 


UMULASSif.IED 


slcurity  classification  or  this  face  rwf»«n  om«  k» 


REPORT  DOCUMENTATION  PAGE 

KCAD  WmtUCTtCMt 

BEFORE  COMPLETING  FORK 

1.  REPORT  NUMBER  1.  GOVT  ACCESSION  NO. 

TECHNICAL  REPORT  ARBRL-TR- 02496  ft  i  ^X 

>.  RECIPIENT'S  CATALOG  NUMBER 

4-  TITLE  <r iblffltj 

A  SURVEY  OF  PENETRATION  MECHANICS  FOR  LONG  RODS 

1.  TYPE  OF  REPORT  A  PERIOD  COVERED 

A.  3ERFORMINO  ORO.  REPORT  NUMBER 

7.  AUTHOR!  •) 

T.  W.  Wright 

B.  CONtRACT  ON  QRANT  NUMBER**) 

•  PERFORMIN'!  ORGANIZATION  NAME  AND  AOORESS 

US  Army  Ballistic  Research  Laboratory 

ATTN:  DRDAR-BLT 

Aberdeen  Proving  Ground,  MD  21005 

10.  PROGRAM  ELEMENT.  PROJECT  TASK 

n.  controlling  office  name  ano  aodress 

US  Army  Armament  Research  §  Development  Command 

US  Army  Ballistic  Research  Laboratory  (DRDAR-BLA-S) 
Aberdeen  Proving  Ground,  MD  21005 

12.  REPORT  DATE 

June  1983 

12.  NUMBER  or  PAf.ES 

43 

14  MONITORING  AGENCY  name  A  AOORCSSflf  dlttmrmnt  hom  Controlling  OtHca) 

IS.  SECURITY  CLASS,  (vl  (Ala  npart) 

UNCLASSIFIED 

tSa.  DECLASSIFICATION/ DORN  OP  ADI  NO 
SCHEDULE 

i«.  Distribution  statement  (oi  t him  Raportj 


Approved  for  public  release;  distribution  unlimited. 


17  OiST  RIBUTION  STATEMENT  ot  (A*  •battmct  m>  taraW  In  Block  W.  If  dlfloronl  tram  Report) 


it.  supplementary  notes 


19.  KEY  WORDS  fContlnu*  on  r«r«r««  «/rf*  It  nocmuimry  mid  identity  by  bloeJt  nuaib«0 

Penetration  mechanics,  long  rods,  armor,  kinetic  energy  penetrators 


*0  ABSTRACT  fGcottaua  an  onto  *MM  H  immi'  tmd  It'wullr  *T  Mock  nypoocj 

Some  of  the  simpler  methods  in  current  use  for  analyzing  ballistic  impacts  by 
long  rod  penetrators  are  reviewed  and  critiqued.  In  spite  of  several  obvious 
shortcomings  the  eroding  rod  model  has  been  used  as  a  starting  point  and  guide 
for  an  experimental  investigation  of  penetration  phenomena.  Some  of  the 
principal  results  of  that  program  are  reviewed,  and  several  areas  for  future 
work  are  identified.  Carefully  developed,  time  resolved  data  of  the  types 
described  in  this  report  should  serve  as  a  standard  of  comparison  for  large- 


DO  ,  jam  71  1^*3  EOtTlOW  or  »  MOV  ••  IS  OBSOLETE  UNCLASSIFIED 


S.  CUBITY  CLMSIFICA'IOM  or  T»ll«  PAGE  Data  CnfaraO 


SMMWTV  CLASSIFICATION  IF  T»IS  FAOSfl»«i 


scale  computations  of  penetration  phenomena.  These  data  are  already  being 
used  to  guide  development  of  improved  engineering  models  for  long  rod 
penetration. 


j 


StCUNITV  CLASSIFICATION  OF  THIS  F  AOtfWlwl  Dm*  SntorMt) 


TABLE  OF  CONTENTS 


LIST  OF  ILLUSTRATIONS . 5 

I.  INTRODUCTION . 7 

II.  DATA  CORRELATION  AND  DIMENSIONAL  ANALYSIS . 7 

III.  ENGINEERING  ANALYSIS  .  13 

IV.  AN  EXPERIMENTAL  PROGRAM . 20 

A.  Experiments  with  Instrumented  Rods . 21 

B.  Experiments  to  Measure  Penetration  .  25 

C.  Experiments  to  Measure  Target  Stress  .  29 

V.  OTHER  CONSIDERATIONS  .  31 

VI.  COMPUTER  CODES . 32 

VII.  CONCLUSIONS . 33 

REFERENCES . 34 


DISTRIBUTION  LIST 


37 


FIG.  NO. 
1 
2 

3 

4 

5 

6 

7 

8 

9 

10 


LIST  OF  ILLUSTRATIONS 


Page 


Principal  features  of  long  rod  penetration .  8 

Limit  velocity  as  calculated  from  Cqn.  (1)  compared  with 

observed  values .  12 

Failure  of  extrapolation  with  a  power  law  .  14 

Typical  predictions  of  the  eroding  rod  model .  16 

Experimental  curve  of  penetration  vs.  striking  velocity 

(after  Stilp  and  Hohler11) . 18 

Curves  of  constant  strain  in  rod  coordinates  for  a  high- 

strength  steel  rod  and  an  armor  steel  target.  Impact  speed 

was  1000  m/s .  22 

Calculated  stress  and  particle  velocity  for  rod  in  Fig.  6  at 
two  impact  speeds.  Also  shown  is  the  measured  quasi-static 
stress/strain  curve  for  the  rod . 24 

Measured  penetration  vs.  time  for  same  rod  and  target  as 

in  Fig.  6 .  26 


Composite  penetration  diagram  in  laboratory  coordinates  ....  28 

Measured  target  stress  vs.  time  at  depths  of  12.7  mm  and 
19.0  mm  below  the  impact  surface.  The  rod  is  a  tungsten 
alloy  and  the  target  is  armor  steel  . 


30 


I.  INTRODUCTION 


Long  rod  penetrators  are  modern  equivalents  of  the  cannonball,  intended 
to  pierce  a  target  by  depositing  large  amounts  of  kinetic  energy  in  a  con¬ 
centrated  region.  Typical  striking  velocities  lie  in  the  range  of  approxi¬ 
mately  1-3  km/s.  Although  a  modern  penetrator  is  a  machined  metal  part  of 
some  complexity,  since  it  must  be  fitted  with  various  appurtenances  to  achieve 
efficient  launch  and  flight,  it  may  be  idealized  as  a  right  circular  cylinder 
with  a  ratio  of  original  length  to  diameter,  L0/D,  greater  than  some  arbitrary 
number,  say  ten.  The  penetrator  is  intended  to  fly  and  to  strike  its  target 
end  on.  When  the  target  thickness  is  greater  than  a  few  penetrator  diameters, 
the  mechanism  of  penetration  is  a  complex  process  in  which  a  cavity  is  made  in 
the  target  and  the  end  of  the  penetrator  erodes  away.  Some  of  the  principal 
features  of  impact  are  shown  in  Figure  1. 

Material  response  properties  under  impact  conditions  are  also  complex. 

It  has  been  estimated  in  a  recent  report  of  the  National  Materials  Advisory 
Board1  that  typical  maximum  pressures  and  maximum  strain  rates  during  penetra¬ 
tion  may  lie  in  the  ranges  5-50  GPa  and  10  -106  s  ,  respectively.  The  higher 
values  are  associated  with  the  earliest  times  after  impact  and  with  shock 
waves,  but  pressures  and  strain  rates,  aided  by  the  presence  of  nearby  free 
surfaces,  quickly  decay  toward  the  lower  values  for  the  bulk  of  the  process. 
Since  strong  armor  or  penetrator  materials  may  have  ultimate  stresses  of 
1-2  GPa,  it  is  apparent  that  maximum  pressures  are  not  orders  of  magnitude 
higher  than  material  strengths  except  perhaps  locally  for  very  short  periods 
of  time.  In  fact,  because  of  wave  reflections,  tensile  stresses  occur  as  well. 
In  short,  although  stresses  are  high,  they  are  not  so  high  that  material  prop¬ 
erties  can  be  safely  ignored.  Furthermore,  phase  changes  may  occur  in  some 
materials,  notably  Fe  and  U,  high  pressures  are  well  known  to  have  a  strong 
influence  on  flow  and  fracture  characteristics,  and  strain  rates  are  high 
enough  to  affect  the  response  of  many  materials. 

On  the  other  hand,  it  is  by  no  means  clear  how  many  details  of  physical 
and  material  processes  are  really  required  in  order  to  achieve  accurate  and 
useful  ballistic  predictions. 

There  are  three  levels  of  analysis  currently  used  in  an  effort  to  come  to 
grips  with  the  difficulties  outlined  above.  In  order  of  increasing  complexity 
in  the  representation  of  physical  and  material  response  phenomena  these  are: 

(1)  empirical  data  correlation,  (2)  engineering  analysis,  and  (3)  large-scale 
computer  codes.  Each  approach  is  useful,  each  has  its  own  limitations,  and  to 
be  effective  each  must  be  supported  by  its  own  brand  of  high  quality  experi¬ 
mentation.  In  this  survey  some  comments  will  be  made  on  all  three  approaches, 
but  emphasis  will  be  placed  on  one  particular  engineering  model  and  related 
experimental  results. 


II.  DATA  CORRELATION  AND  DIMENSIONAL  ANALYSIS 

At  the  simplest  engineering  level  the  ballistician  would  like  to  be  able 
to  predict  the  outcome  of  an  encounter  between  a  penetrator  and  target. 


1  "Vtateviale  Response  to  Ultra-High  Loading  Rates,"  Rational  Materials  Advisory 
Board  Publ.  RMAB-8S6,  Washington:  Rational  Academy  of  Sciences, 1980. 
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Figure  1.  Principal  features  of  long  rod  penetration. 


That  is,  with  given  information  about  geometries  and  material  identifications 
of  both  penetrator  and  target  and  on  line  of  attack  and  striking  velocity,  one 
would  like  as  a  minimum  to  predict  whether  or  not  the  target  will  be  perfo¬ 
rated,  one  would  like  to  calculate  depth  of  penetration,  and  if  it  is  perfo¬ 
rated,  to  calculate  residual  velocity  and  residual  mass.  Consequently,  ballis¬ 
tic  testing  today  concentrates  on  careful  control  of  line  of  flight  and  target 
attitude,  and  on  accurate  measurement  of  striking  and  residual  velocities. 

When  flash  X-rays  are  used  for  velocity  measurements,2  information  is  also 
gained  on  projectile  yaw  at  impact  and  line  of  flight  after  impact.  Often  an 
estimate  of  residual  mass  can  be  made  as  well. 

There  have  been  many  attempts  to  make  empirical  correlations  of  ballistic 
data.3  Typical  cases  are  listed  below  in  Table  I. 


TABLE  I:  Typical  Correlation  Functions  for  Ballistic  Data 
deMarre4  :  MV^2  *=  CDfi  f* 

Grabarek5:  MV^S  *  C 

Thor6  :  Vg  -  VR  «  C  T®  VgY 

M  «  projectile  mass  Vg  ■  Striking  velocity 

D  »  projectile  diameter  VR  ■  residual  velocity 

T  »  target  thickness  VL  »  limit  velocity 

Q  »  obliquity,  measured  from  plate  normal 
a,  B,  Yi  C  »  empirical,  best  fit  parameters 


V.  Grabarek  and  L.  Herr,  "X-Ray  Multi-Flash  System  for  Measurement  of  Pro¬ 
jectile  Performance  at  the  Target ,"  BRL-TN-1634,  Ballistio  Research  Laboratory 
APG,  MD,  1966  (AD  807619). 

3M.  E.  Bachman  and  W.  Goldsmith,  "! The  Mechanics  of  Penetration  of  Projectiles 
into  Targets,"  Int.  J.  Eng.  Sci. ,  16,  1-99,  1978. 

*V.  deMarre,  "Perforation  of  Iron  and  Steel  Sheets  with  Normal  Firing  (Trans.), 
Memorial  de  I’Artillerie  de  Marin 4  14,  1886.  Empirical  formulas  of  this  type 
for  limit  energy  are  often  called  deMarre  relations ,  at  least  in  the  V.S.  In 
his  original  paper  deMarre  gave  the  formula  with  a  *  1.4  and  B  m  l.S. 

5 <7.  Grabarek,  "Penetration  of  Armor  by  Steel  and  High  Density  Penetrators 
BRL-MR-2134,  Ballistic  Research  Laboratory ,  APG  MD,  1971  (AD  618394L). 

6  Project  Thor  Tech.  Report  No.  47,  " The  Resistance  of  Various  Metallic  Mater¬ 
ials  to  Perforation  by  Steel  Fragments ;  Empirical  Relationships  for  Fragment 
Residual  Velocity  and  Residual  Weight,"  Ballistic  Analysis  Laboratory,  The 
Johns  Hopkins  University,  1961. 
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These  correlations  and  others  like  then  can  be  extrwely  useful  for  pre¬ 
dicting  other  cases  that  lie  within  the  same  general  data  set,  that  is,  for 
data  interpolations  or  even  for  small  extrapolations.  It  will  be  noted,  how¬ 
ever,  that  none  of  the  three  cases  cited  contains  information  on  material 
properties.  Therefore,  every  change  in  material  roquires  a  new  set  of 
coefficients.  In  fact,  none  of  the  formulae  in  Table  I  is  even  given  in  non- 
dimensional  form,  so  that  the  range  of  applicability  in  each  case  is  unduly 
restrictive. 

On  the  other  hand,  there  are  many  dependent  variables  of  possible  interest, 
and  even  more  independent  variables  that  may  influence  the  outcome  of  a  given 
encounter.  Some  of  the  possibilities  are  shown  in  Table  II.  The  rows  of  dots 
signify  all  other  possible  variables,  and  each  entry  in  the  list  of  material 
variables  can  be  given  for  both  penetrator  and  target.  Even  with  the  reduction 
to  nondimensional  variables,  the  analyst  is  still  faced  with  a  bewildering 
array  cf  possible  influences,  and  because  of  limitations  in  both  time  and  money 
available  for  testing,  it  invariably  turns  out  that  the  test  data  are  incom¬ 
plete  with  respect  to  some  variables . 


TABLE  II:  Variables  for  Dimensional  Analysis 


Dependent  Variables: 

P  *  depth  of  penetration 
VL  ■  limit  velocity 
VR  *  residual  velocity 
Lp  »  residual  length 

Independent  Variables: 

Physical  Characteristics 

Lq  =  initial  length 
D  =  diameter 

Kinematic  Characteristics 

Vg  =  striking  velocity 
o  =  angle  of  yaw 

Material  Characteristics 

Ey  *  yield  stress 

E  *  ultimate  stress 
u 

h  *  hardness 
■  elongation  at 


MR  -  residual  mass 

9R  *  angle  of  residual  trajectory 

vh  »  hole  volume 


=  initial  mass 
o 

T  =■  target  thickness 


0 g  *  angle  of  striking  trajectory 


p  ■  density 

Kjc  ■  fracture  toughness 
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Recently  an  attempt  has  been  made  by  Bruchey7  to  include  more  variables 
than  has  been  customary  in  the  past.  The  nondimensional  form  chosen  by  him 
for  data  correlation  is  given  in  F.qn.  (1) . 


E  and  E_  both  have  the  form 
P  t 

E  “  1  (ZY  +  V  (2) 

and  represent  measures  of  the  energy  per  unit  volume  that  can  bo  absorbed  in 
the  penetrator  and  target  materials  before  rupture  occurs.  Tsece  is  the  line 
of  sight  thickness  of  the  target,  and  all  other  quantities  are  as  listed  in 
Table  II.  The  coefficient  C  and  the  powers  a,  B,  and  y  are  determined  by  a 
multilinear  regression  analysis.  The  total  data  set  includes  only  25  separate 
cases,  but  within  that  set  there  is  substantial  variation  of  the  nonlimensional 
ratios.  The  ratio  Lq/D  varies  by  a  factor  of  2,  TsecG/D  varies  by  4,  Pp/Pt 

varies  by  1  1/4,  E^/E^  varies  by  5,  and  varies  by  4. 

Other  groupings  of  the  variables  could  have  been  made.  For  example, 
instead  of  the  left-hand  side  of  Eqn.  (1)  as  given,  the  grouping 


E^T  secG 

could  have  been  used.  This  form  may  be  interpreted  as  the  kinetic  energy 
available  per  unit  area  normal  to  the  line  of  flight  compared  with  the  ability 
of  the  target  to  absorb  energy  per  unit  area.  When  this  form  is  used,  the 
powers  on  the  right  hand  side  are  all  small  numbers  less  than  one,  indicating 
only  weak  dependence,  but  the  overall  agreement  after  the  regression  analysis 
is  no  better  than  the  original  result. 

For  Eqn.  (1)  the  calculated  vs.  the  observed  limit  velocities  are  shown 
in  Fig.  2. 

For  velocity,  varying  by  a  factor  of  2,  agreement  is  generally  within 
about  10%.  For  more  restricted  data  sets,  that  is  for  data  sets  with  less 
variation  in  material  and  geometric  properties,  much  better  agreement  can  be 
achieved,  but  only  with  a  corresponding  loss  in  generality  of  the  correlations. 

In  that  it  attempts  to  incorporate  seme  effect  of  material  properties, 

Eqn.  (1)  is  an  improvement  over  Table  I,  but  it  still  uses  only  a  few  of  the 
non-dimensional  numbers  possible  from  Table  II.  Furthermore,  the  assumed 


W.  Bruchey ,  Private  communication.  In  preparation  ae  a  report ,  Ballistic 
Research  Laboratory 4 PC,  MD. 
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power  law  relationship  among  variables  is  an  extremely  limited  choice  of 
possible  functional  dependences. 


As  illustrated  in  Fig.  3  extrapolation  with  the  weak  power  law  dependence 
may,  and  in  fact  often  does,  cause  rubstantial  error.  On  the  other  hand, 
experimentation  to  determine  genera*  ,'ependence  of  the  form 


2 

ppVL  Lo  f  / T  sec0  Ep  Pp  Eo 
E  T  sec©  \  D  *  E  ’  p  .  *  D 


(3) 


for  example,  would  be  an  enormous  task.  A  better  strategy  would  be  to  look 
for  trends  within  subsets  of  the  available  data,  where  only  one  variable  at  a 
time  is  allowed  to  vary,  and  then  to  try  to  fit  the  data  with  simple  function 
that  have  limited  degrees  of  freedom. 


But  no  matter  how  the  correlations  are  made,  it  will  be  difficult,  if  not 
impossible,  to  gain  more  than  qualitative  insight  into  the  nature  of  the 
physical  processes  at  work.  For  a  quantitative  view  of  these  processes,  it  is 
necessary  to  construct  physical  models. 


III.  ENGINEERING  ANALYSIS 

There  have  been  many  attempts  to  construct  simple  models  of  penetration. 
Many  of  these  are  described  in  a  recent  survey,3  but  only  one,  the  eroding  rod 
model,  will  be  summarized  here,  because  it  has  proved  useful  as  a  starting 
point  for  an  experimental  program  and  for  further  theoretical  refinements. 

The  goal  for  any  engineering  model  is  to  construct  a  system  of  equations  with 
enough  detail  to  be  physically  realistic  but  at  the  same  time  simple  enough 
for  rapid  computation. 

The  eroding  rod  model  in  its  simplest  form  was  given  independently  by 
Alekseevskii8  and  Tate.9  With  the  instantaneous  length  of  the  penetrator  L, 
speed  U,  and  depth  of  penetration  P,  the  equations  may  be  written  as  follows: 


Mass: 

L  =  P  -  U 

(4) 

Momentum : 

Pp  LU  =  -Zp 

(5) 

Modified  Bernoulli: 

ZP  +  h  ep(u-p)2  «  +  i  pt  *2 

(6) 

In  these  equations  the  only  material  properties  included  are  flow  stress  E  and 
density  P  where  subscripts  p  and  t  denote  penetrator  and  target,  respectively. 


eF.  P.  Alekseevskii.  "Penetration  of  a  Bod  into  a  Target  at  High  Velocity ," 
Comb .,  Expl.,  and.  7nock  Waves,2,  63-66 ,  1966 ,  trans.  from  Russian. 

9A.  Tate,  "A  Theory  for  the  Deceleration  of  Long  Rods  After  Impact , "  J.  Meah. 
Phys.  Sol.  15,  387-399,  1967,  and  "Further  Results  in  the  Theory  of  Long  Rod 
Penetration,"  J.  Meah.  Phys.  Sol.,  17,  141-150,  1969. 
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The  interpretation  of  Equations  (4)  and  (5)  as  balance  of  mass  and  balance 
of  linear  momentum  is  straight  forward.  Equation  (6)  should  be  interpreted  as 
shewing  equality  of  axial  stresses  at  the  interface  between  penetrator  and 
target  on  the  stagnation  line.  For  the  case  of  steady  flow  with  perfect  fluids 

1  2 

the  expression  p  +  j  pv  is  constant  along  a  steamlire.  At  the  stagnation 

point  the  velocity  relative  to  that  point  vanishes  so  the  constant  is  the 
interface  pressure.  If  both  the  rod  and  target  are  assumed  to  behave  like 
perfect  fluids  in  steady  flow  for  some  short  distance  away  from  the  interface 
where  the  solid  character  of  each  material  reasserts  itself,  then  the  left 
side  of  (6)  holds  for  the  rod  and  the  right  side  for  the  target  with  the  equal¬ 
ity  holding  because  stresses  must  match  at  the  interface. 

Inherent  in  the  model  is  the  assumption  that  both  penetrator  and  target 
may  be  approximated  as  rigid/plastic  materials.  That  is,  either  material 
flows  at  its  characteristic  stress  E  or  it  locks  up  and  doesn't  deform  at  all. 
For  example,  if  lt  >  Ep,  both  penetrator  and  target  material  flow  only  if 

U  >  ^2(Et  ~  ^p^  /pp  ‘  At  this  critical  speed  P  becomes  zero  and  the  target 

ceases  to  deform,  but  the  penetrator  continues  to  erode  until  U  drops  to  zero. 
On  the  other  hand,  if  £p  >  Et ,  both  materials  flow  only  if  P  >^2(Ep  -  l^)/  pt* 

When  the  rate  of  penetration  drops  below  that  speed,  erosion  of  the  penetrator 

•  1  *2 
stops,  U  =  P,  and  the  right  side  of  (5)  must  be  replaced  by  -(E^  +  ^f^P  ) • 

The  penetrator  then  continues  to  penetrate  as  a  rigid  body  until  U  drops  to 
zero . 


Equations  (4)  -  (6)  are  simple  enough  to  be  solved  on  a  prograranable  hand 
calculator.  Typical  predictions  are  shown  in  Fig.  4  where  nondimens ional 
penetration  is  plotted  as  a  function  of  nondimensional  striking  velocity. 

The  qualitative  features  of  the  model  are  striking  and  useful  in  themselves. 
The  curve  moves  to  the  right  for  increasing  target  strength  and  to  the  left 
for  increasing  penetrator  strength,  which  is  in  accord  with  common  sense,  of 
course,  but  the  model  gives  a  quantitative  estimate  of  the  amount  of  change 
as  well.  All  curves  have  a  characteristic  S-shape  and  show  a  saturation  value 
for  high  impact  speeds.  This  is  in  sharp  contrast  to  the  shape  of  the  power 
law  curves  customarily  used  for  data  correlation.  The  eroding  rod  model  also 
brings  out  important  natural  scaling  effects,  the  most  important  one  being 
that  penetration  is  proportional  to  penetrator  length  if  all  other  properties 
are  held  fixed.  Thus,  to  maximize  penetration  for  a  penetrator  of  fixed  mass, 
the  ratio  LQ/D  should  be  as  high  as  possible,  consistent  with  structural 

stability  and  other  engineering  constraints. 

A  less  obvious  consequence  of  the  model  is  that  for  a  fixed  amount  of 
kinetic  energy  and  fixed  LQ/D  there  is  an  optimum  striking  velocity  that  will 

maximize  penetration.  That  same  velocity  also  minimizes  the  energy  required 
to  achieve  a  fixed  depth  of  penetration.  These  facts  follow  simply  from  the 
saturation  effect.  Penetration  may  be  represented  as  P  *  Lq  S(Uq)  where 

S(Uq)  is  the  S-shaped  saturation  curve.  Then,  since  kinetic  energy  K  is 
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3  2 

proportional  to  lQ  Uq  for  fixed  1,Q/D,  the  results  follow  from  a  minimization 
problem  with  side  constraint.  The  optimum  velocity  is  the  solution  of 


dS  2  S_ 

<fiT  ‘  3  U  ’ 
o  o 


(7) 


and  may  easily  be  found  graphically.  There  will  be  no  solution  to  (7)  if  the 
curve  S(U  1  does  not  bend  over  far  enough  so  the  tendency  toward  saturation  is 

essential  for  this  result  to  hold.  Results  of  this  kind  have  previously  been 
given  by  Frank10 

Figure  5  shows  experimental  data  reported  by  Stilp  and  Hohler11  Similar 
data  were  given  oy  Tate,  et  at  12  The  striking  simularity  of  the  experimental 

curve  and  the  theoretical  curve  for  l  <  E*  is  clearly  evident. 

P  t 

Equations  (4)  -  (6),  therefore,  have  considerable  intuitive  appeal,  com¬ 
bining  as  they  do  simplicity  with  qualitative  accuracy,  but  they  are  difficult 
to  use  quantitatively  because  values  for  the  characteristic  flow  stresses  are 
not  readily  available  a  priori.  The  usual  procedure  has  been  to  choose  the 
stresses  so  as  to  fit  an  experimental  S-curve,  and  then  to  consider  them  as 
material  constants.  Due  to  approximations  inherent  in  the  model,  this  pro¬ 
cedure  is  unlikely  to  be  satisfactory  in  all  cases.  There  are  several  diffi¬ 
culties,  nearly  all  of  which  are  associated  with  the  modified  Bernoulli  equa¬ 
tion  itself.  These  have  been  described  in  seme  detail  previously  by  Wright13 

The  first  problem  is  in  the  origin  of  Equation  (6).  Strictly  speaking, 
the  equation  only  applies  for  steady  flow  in  perfect  fluids.  The  stresses  in 
the  equation  would  be  those  that  occur  at  interfaces  where  the  penetrator  or 
target  material  undergoes  a  transformation  from  rigid  solid  to  perfect  fluid. 
Furthermore,  the  fluid  layers  must  be  imagined  to  have-  constant  thickness  so 
that  the  fluid/solid  interfaces  both  have  the  speed  P.  Melting  can  occur  in 
ballistic  impact  if  the  speed  is  high  enough,  but  in  any  case  there  will  be  a 
region  of  plastic  flow  on  either  side  of  the  penetrator/target  interface.  For 
normal  impact  with  the  z-axis  along  the  centerline  of  the  rod,  the  z-component 
of  the  momentum  equation  may  be  written 


10K.  1 Frank,  "A  Qualitative  Determination  of  the  Velocity ,  Mass,  and  Fineness 
Ratio  Required  to  Defeat  Single  Plate  Targets ,"  Spring  Tech.  Conf . ,  Ballistic 
Research  Laboratory,  APG,  MD,  1978. 

11  7.  Hohler  and  A.  J.  Stilp,  " Penetration  of  Steel  and  High  Density  Rods  in 
Semi-Infinite  Steel  Targets,  "  Third  Int.  Sijmp.  Ballistics,  Karlsruhe, 

Germany,  1977. 

124.  Tate,  K.  E .  B.  Green,  P.  G.  Chamberlain,  and  R.  G.  Baker,  " Model  Scale 
Experiments  on  Long  Rod  Penetrators , "  Fourth  Int.  Synp.  Ballistics,  Monterey, 
CA,  1978. 

13 7.  W.  Wright,  "Penetration  with  Long  Rods:  A  Theoretical  Framework  and  Com¬ 
parison  wi.th  Instrumented  Impacts,  "  ARBRL-TR-02323,  Ballistic  Research 
Laboratory,  APG,  MD,  1981  (AD  A101344). 
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where  the  commas  denote  partial  differentiation.  On  the  centerline  the  trans¬ 
verse  components  of  velocity,  u  and  u  ,  are  zero  by  symmetry  and  for  steady 

x  y 

flow  is  zero,  also.  Shear  stresses  also  vanish  on  the  centerline,  but 

Z  1 1 

not  their  derivatives,  so  integration  of  Equation  (8)  between  points  a  and  b, 
assuming  incompressibility,  yields 


t 

zz 


(b)  -  t22(a) 


*  Sx.x5  dz 


(9) 


If  point  b  is  the  stagnation  point  at  the  penetrator/target  interface,  and  if 

point  a  lies  at  the  rigid/plastic  boundary  in  either  penetrator  or  target, 

then  with  t  (a)  =  -I  for  compressive  stresses  Equation  (6)  must  be  replaced  by 
z  z 

V  ViVu4)2  -Et  *  h  cioj 

wher-  I  is  the  integral  ir.  Equation  (9).  If  l  and  I  are  lumped  together,  then 
Equation  (10)  has  the  same  form  as  Equation  (6),  but  the  integrals  are  undeter¬ 
mined  . 


A  second  problem  concerns  the  validity  of  the  assumptions  of  steady  flow 
and  rigid/plastic  materials.  Taken  overall,  kinetic  energy  penetration  is 
clearly  an  unsteady  process,  espec-’lly  in  the  initiation  and  final  stopping 
or  breakout  phases,  but  there  may  be  parts  of  the  process,  localized  in  both 
space  and  time  that  are  nearly  steady.  For  example,  during  the  intermediate 
stages  of  penetration  a  local  interaction  region  may  form  near  the  penetrator/ 
target  interface,  within  which  Equation  (10)  holds  approximately.  If  the 
target  is  thick  enough  and  the  penetrator  long  enough,  this  intermediate  stage 
will  be  the  dominant  phase  of  penetration. 

Even  in  such  a  case,  however,  there  is  no  clear  way  to  make  an  accurate 
estimate  of  the  characteristic  stresses  from  knowledge  of  measured  mechanical 
properties.  The  stress/strain  curve  of  most  real  high-strength  materials 
cannot  be  well  approximated  by  the  rigid/plastic  assumption  because  of  work 
hardening.  Furthermore,  failure  processes  in  a  high  rate  interaction  zone 
must  be  regarded  as  largely  unknown  so  that  it  is  impossible  to  specify  either 
a  maximum  strain  or  a  maximum  stress  at  failure  with  any  degree  of  confidence. 
Thus,  there  is  at  present  no  well  defined  way  to  choose  the  characteristic 
flow  stresses.  Even  a  criterion  based  on  energy  equivalence  does  not  seem 
particularly  feasible. 

In  addition  to  the  uncertainty  concerning  flow  stresses  and  the  integrals 
in  Equation  (10) ,  there  is  further  uncertainty  concerning  the  spherical  com¬ 
ponent  of  stress  in  both  penetrator  and  target.  The  axial  stress  in  Equation 
(9)  or  (10)  may  be  decomposed  into  the  sum  of  a  spherical  and  a  deviatoric 
part.  The  flow  stress  is  associated  only  with  the  latter,  but  the  spherical 
or  hydrostatic  part  is  totally  undetermined.  In  the  penetrator  nearby  free 
lateral  surfaces  will  tend  to  minimize  the  spherical  stress,  but  in  the  target 
there  are  no  lateral  surfaces,  and  the  hydrostatic  component  will  be  the  major 
part  of  the  axial  stress. 
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Still  another  consequence  of  real  nonrigid  Material  behavior  is  that 
there  will  be  considerable  plastic  deformation  at  some  distance  from  the  quasi¬ 
steady  Interaction  zone.  In  the  rod  this  has  the  effect  of  slowing  material 
down  before  it  reaches  that  zone  so  that  the  speed  U  to  be  used  in  Equation  (10) 
is  less  than  that  of  the  rear  end.  In  the  target,  plastic  deformation  (includ¬ 
ing  bulging  of  the  rear  surface)  has  the  effect  of  increasing  the  spatial  rate 
of  penetration.  Target  deformation  and  bulging  has  been  explicitly  included 
in  the  model  for  perforation  by  rigid  projectiles  given  by  Ravid  and  Bodnev.1* 

Finally,  inertia  will  tend  to  increase  the  crater  size  even  after  the 
active  driving  forces  have  ceased.  The  effect  of  crater  inertia  is  particu¬ 
larly  pronounced  for  low  ratios  of  L^/D  or  ultra  high  velocities  and  has  the 
affect  of  increasing  penetration  above  the  curves  calculated  in  Fig.  4.  This 
effect  has  been  considered  by  Frank15  and  more  recently  by  Tate.16 

The  third  and  final  problem  is  that  the  stress  tending  to  decelerate  the 
rod  in  Equation  (5)  is  an  average  stress  over  the  cross  section,  whereas  the 
stress  that  enters  into  the  modified  Bernoulli  equation  is  the  local  stress  on 
the  centerline.  These  may  be  somewhat  different  from  each  other  especially  if 
the  integral  in  Equation  (9)  is  lumped  together  with  the  axial  stress  to  give 
an  effective  value  for  £. 

For  all  of  the  above  reasons  -  the  origin  and  meaning  of  the  modified 
Bernoulli  equation,  the  questionable  validity  of  steady  flow  and  rigid/plastlc 
assumptions,  and  the  inequality  of  average  and  centerline  stresses  -  the 
eroding  rod  model  is  difficult  to  use  quantitatively  and  gives  only  limited 
insight  into  the  details  of  the  actual  penetrator/target  interaction.  Never¬ 
theless,  in  spite  of  its  shortcomings  the  model  appears  to  contain  the  germ  of 
a  sound  theory  of  penetration  for  long  rods,  as  shown  by  the  qualitative 
success  of  its  predictions. 


IV.  AN  EXPERIMENTAL  PROGRAM 

Penetration  is  controlled  by  flow  and  failure  processes  that  occur  at  or 
near  the  penetrator/target  interface.  These  processes  depend  only  on  the 
material  properties  of  the  particular  materials  involved,  and  therefore,  in 
large  measure  they  will  depend  only  on  the  local  stress  and  strain  fields,  and 
possibly  the  gradients  or  time  rates  of  change  of  those  fields.  These  ideas 
are  at  least  compatible  with  the  formation  of  a  region  of  quasi-steady  flow  at 
the  interface.  Since  it  is  a  local  phenomenon,  it  seems  reasonable  to  suppose 
that  the  steady  conditions  established  will  depend  only  on  the  local  geometry 
and  the  mass  flux  through  the  region. 


l4W.  Ravid  and  S.  R.  Bodner,  "Dynamic  Perforation  of  Vie  coplastic  Plates  by 
Rigid  Projectiles , "  Report  of  the  Material  Mechanics  Laboratory ,  Teohnion, 
Haifa,  Israel,  1982. 

1SK.  Frank,  Unpublished  Vork.  His  essential  idea  was  to  combine  the  results 
of  the  eroding  rod  model  with  hypervelocity  cratering  results  of  Christman 
and  Gehring,  J.  Appl.  Phy s.,  37,  1579-1587,  1966. 

1 6  A.  Tate,  " Extensions  to  the  Modified  Hydrodynamic  Theory  zf  Penetration, " 
Preprint,  to  be  published.  20 


With  the  ideas  of  the  preceding  paragraph  as  a  working  hypothesis  several 
experiments  have  been  devised  to  begin  probing  the  details  of  the  interaction 
process.  Obviously,  measurements  cannot  be  made  directly  in  the  interaction 
region,  but  they  can  be  made  in  both  penetrator  and  target  in  material  immedi¬ 
ately  adjacent  to  it.  In  fact,  in  terms  of  the  modified  Bernoulli  equation, 
it  has  been  possible  to  make  independent  measurement  of  E^,  Et>  U,  and 

A.  Experiments  with  Instrumented  Rods. 

In  the  first  of  these  experiments,  done  by  Hauver17*18  and  Hauver  and 
Melani,  strain  gages  are  placed  on  the  long  rod  penetrators  and  strain-time 
histories  are  recorded  at  several  stations  on  the  rod.  Then  from  the  one¬ 
dimensional  theory  of  wave  propagation  it  is  possible  to  reconstruct  the  stress 
and  particle  velocity  histories  throughout  the  rod.  This  is  done  as  follows. 
The  simplest  equations  for  1-D  wave  propagation  may  be  written 

S»z  »  pw,t  ;  w,2  -  e,t  (11) 

where  E  is  the  average  axial  stress  in  the  rod  (referred  to  the  initial  cross 
section),  p  is  the  initial  density,  W  is  the  particle  velocity,  and  e  is 
engineering  strain.  The  spatial  location,  z,  is  a  function  of  the  material 
particle,  Z,  and  time,  t,  so  that 

z  »  z(Z.t),  N  *  2,t  e  »  x»2  -1  (12) 

With  e-t  data  for  several  stations  at  hand  curves  of  constant  strain  may  be 
plotted  in  Z-t  space.  Along  such  a  curve  the  slope,  c  ,  called  the  plastic 
wave  speed,  may  be  measured.  p 

de  «  e,  dt  +  c >zdZ  ■  0  (13) 


dZ 

dt 


(14) 


These  curves  need  not  be  straight  lines,  but  if  they  are,  as  all  experiments 
to  date  indicate,  the  analysis  is  greatly  simplified  because  c^  then  depends 

only  on  e.  Typical  data  is  shown  in  Fig.  6,  taken  from  Hauver.18  The  actual 
experiments  were  performed  in  a  4”  light  gas  gun  with  the  penetrator  stationary 
for  ease  of  instrumentation  and  the  target  plate  launched  from  the  4"  gun. 

The  experiment  is  shown  schematically  in  the  figure. 


17C.  E .  Hauver,  " Penetration  with  Instrumented  Rods,"  Int.  J.  Eng.  Soi. ,  16, 
871-87?,  1978. 

1 8  G.  E.  Hauver,  "Experiments  with  Instrumented  Long-Rod  Penetrators,"  Fifth 
Int.  Symp.  Ballistics,  Toulouse ,  France,  1980. 

i9G.  £’.  Hauver  and  A.  Melani,  "Strain-Gage  Techniques  for  Studies  of  Projectiles 
During  Penetration ,"  ARBRL-MR-03082 ,  Ballistic  Research  Laboratory,  APG,  MD, 
1981  (AD  A098660). 
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Figure  6.  Curves  of  constant  strain  in  rod  coordinates  for  a  high-strength 
steel  rod  and  an  armor  steel  target.  Impact  speed  was  1000  m/s. 
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Equations  (11)  may  be  intograted  with  respect  to  1  along  lines  t  *  con¬ 
stant  as  has  been  described  by  Kolsky20  and  others.  If  the  rod  has  speed  WQ 
and  zero  stress  and  strain  at  any  station  before  a  wave  arrives,  we  have 


/(z.t) 

W(Z,t)  «  Wo  -f  Cp(c)  de,  (15) 


Z(Z,t) 


(Z,t) 


(c)  de  . 


(16) 


Note  that  the  stress  can  be  obtained  from  Equation  (16)  without  recourse  to 
any  constitutive  assumption.  Typical  results  obtained  by  Hauver17 »18  and 
Hauver  and  Melani19  are  shown  in  Fig.  7. 


The  most  striking  and  immediately  obvious  feature  of  these  results  is 
that  the  stress,  as  ccmputed  from  experiment,  does  not  agree  with  the  measured 
static  stress  and  in  fact  the  computed  stress  depends  on  the  striking  velocity, 
V.  On  the  other  hand,  the  particle  velocity,  as  computed  from  the  experiment, 
agrees  well  with  other  experiments  also  conducted  by  Hauver  and  Melani,19  in 
which  particle  motion  was  determined  directly  by  observing  fiducial  lines  on 
the  penetrator  with  a  streak  camera.  It  is  not  possible  to  measure  stresses 
directly,  but  within  the  limitations  of  Equations  (11),  there  is  no  reason  to 
believe  that  computed  stresses  are  less  accurate  than  computed  particle  veloc¬ 
ities  . 


Two  reasons  have  traditionally  been  given  for  discrepancies  between 
static  and  dynamic  determinations  of  stress/strain  relations.  The  first, 
viscoplastic  or  rate  effects,  has  been  ruled  out  as  a  major  contributor-  for 
the  penetrator  material  in  this  case.  Standard  Hopkinson  bar  tests,  as 
reported  by  Chou?1  have  not  shown  significant  rate  effects.  The  second 
possibility,  the  effect  of  finite  radius,  must  therefore  be  explored. 

Motivated  by  the  results  described  above,  Wright22  has  examined  a  higher 
order  rod  theory  where  radial  strain,  u,  is  treated  as  an  internal  variable, 
and  given  equal  status  with  the  axial  displacement,  W,  as  a  dependent  variable. 
Either  by  use  of  a  variational  principle  to  arrive  at  a  one -dimensional  formu¬ 
lation  directly,  or  by  integration  of  the  three-dimensional  equations  through 
a  cross  section,  one  can  arrive  at  a  coupled  pair  of  one-dimensional  equations. 


20  tf.  Kolsky ,  "Stress  Waves  in  Solids ,"  New  York,  Dover  Publ.,  1963. 

21  P. -C.  Chou,  Letter  Report  to  Ballistic  Research  Laboratory  from  Amy 
Mechanics  and  Materials  Research  Center,  1980. 

22  T.  W.  Wright,  ", Nonlinear  Waves  in  Rods,"  Proc.  IUTAM  Symp.  on  Finite 
Elasticity,  D.  E.  Carlson  and  R.  T.  Shield  (eds.).  The  Hague,  Boston , 
London,  Nijhoff  Publishers ,  1990. 
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Figure  7.  Calculated  stress  and  particle  velocity  for  rod  in  Fig.  6 
at  two  impact  speeds.  Also  shown  is  the  measured  quasi - 
static  stress/strain  curve  for  the  rod. 


PARTICLE  VELOCITY  (u),  mm/fjLS 


3S 

32 


3‘W 

p  T*  ' 

(17) 

3t 

1  oa2  32u 

2pa  “I* 

(18) 

3t 


Here  S  is  the  average  axial  stress,  P  is  the  average  lateral  pressure  on 
the  rod,  Q  is  the  average  polar  moment  of  radial  shear  stress,  a  is  the  initial 
radius  of  the  rod,  and  P,  u,  and  W  are  as  previously  defined.  In  the  constitu¬ 
tive  equations  for  S,  P,  and  Q  the  basic  kinematic  variables  are  the  axial 

strain,  ^  »  the  radial  strain,  u,  and  the  gradient  of  radial  strain, 

Analysis  o 2  these  equations  yields  quantitative  and  qualitative  information 
about  the  structure  of  wave-like  solutions,  either  for  the  nonlinear  elastic 
case  as  in  Reference  22  or  for  the  elastic/plastic  case,  which  is  currently 
being  studied. 

Note  that  Equation  (17)  is  the  same  as  Equation  (1 1) ^ .  Furthermore, 
Equation  (11)2  applies  as  a  compatibility  condition  in  the  present  context  as 
well  as  previously.  Therefore,  in  the  treatment  of  experimental  data  Equa¬ 
tions  (15)  and  (16)  are  still  valid  for  the  computation  of  particle  velocity 
and  stress  in  the  wave.  Now,  however,  there  should  be  no  surprise  that  the 
stress,  as  computed  for  the  dynamic  case,  does  not  agree  with  the  static 
curve  since  the  radial  strain  is  still  completely  undetermined.  Analytically 
it  can  only  be  determined  through  the  coupling  with  Equation  (16),  which 
describes  the  radial  motion  of  the  rod. 

Since  stress  is  also  constant  along  the  lines  of  constant  strain  in 
Fig.  6,  it  would  be  possible  to  determine  the  load  history  on  the  end  of  the 
eroding  penetrator  provided  that  extent  of  erosion  is  known  as  a  function  of 
time  and  that  the  erosion  trajectory  cuts  across  the  straight  lines.  No 
experiment  has  yet  been  devised  to  measure  extent  of  erosion  directly,  but  an 
indirect  approach  seems  possible  by  determining  target  penetration  as  a  func¬ 
tion  of  time  and  matching  that  curve  with  the  projected  arrival  trajectories 
of  various  rod  stations. 

B.  Experiments  to  Measure  Penetration. 

Using  the  same  penetrator  and  target  materials  and  the  same  impact  condi¬ 
tions  as  were  used  by  Hauver,  Netherwood23  attempted  in  a  second  series  of 
experiments  to  measure  rate  of  penetration  by  means  of  insulated  switch  wires 
which  were  inserted  in  small  holes  drilled  into  the  side  of  the  target  block. 
The  situation  is  sketched  in  Fig,  8.  The  penetrator  is  launched  by  a  4"  gas 
gun,  and  the  switch  wires  short  out  in  order  upon  arrival  of  the  penetrator 
at  their  target  stations.  Typical  data  are  shown  in  Fig.  8  with  the  upper 
and  lower  curves  representing  bounds  on  the  data  over  several  repetitions  of 
the  test  with  switch  pins  located  at  different  stations  for  different  tests. 
Distance  is  measured  with  respect  to  the  undeformed  target  and,  therefore,  is 


23P.  H.  Netherwood ,  " Rate  of  Penetration  Measurements*"  ARBRL-MR-02978, 
Ballistic  Research  Laboratory ,  APG *  MD *  1979  (AD  A080S41). 
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a  material  coordinate.  The  curves  show  a  rapid  entry  zone,  followed  by  a 
nearly  steady  zone  of  slow  penetration.  The  last  part  of  the  curves  should 
not  be  interpreted  as  a  rapid  exit  zone,  but  rather  as  an  erratic  failure  zone 
where  a  small  plug  forms  near  the  rear  target  surface  and  shears  off  the  switch 
pins  before  actual  arrival  of  the  penetrator.  The  erratic  behavior  is  believed 
to  be  caused  by  the  somewhat  random  nature  of  material  failure. 

The  detailed  shape  of  the  curves  in  Fig.  8  cannot  be  taken  too  literally 
because  plastic  flow  of  target  material  may  well  cause  the  switch  pins  to 
trigger  before  penetrator  arrival.  Clearly,  however,  the  real  curve  must  lie 
on  or  below  the  apparent  curve  so  that  the  data  give  an  upper  bound  for  depth 
of  penetration  at  any  time.  The  straight,  middle  portion  of  the  apparent 
curve  might  be  supposed  to  lie  parallel  to  the  real  curve  where  quasi-steady 
flow  occurs,  but  doubt  has  been  cast  on  that  interpretation,  as  well,  by 
Hauver's  recent  work.18  It  seems  fair  to  say  that  truly  definitive  measure¬ 
ments  of  penetration  vs.  time  have  yet  to  be  made,  although  Netherwood's 
experiments  represent  a  major  step  in  that  direction. 

To  compare  the  penetration  curve  with  the  projected  arrival  trajectories 
of  various  rod  stations  it  is  necessary  to  use  a  common  coordinate  system. 

The  rod  data  were  taken  in  the  material  coordinates  of  the  rod,  and  the  target 
data  were  taken  in  the  material  coordinates  of  the  target,  but  both  sets  of 
data,  as  well  as  the  curves  of  constant  strain,  may  be  transformed  into  a 
common  spatial  or  laboratory  coordinate  system. 

In  material  coordinates  for  the  penetrator,  the  trajectorie  of  constant 
strain  may  be  written  as  Z  =  Z(t)  so  from  use  of  Equation  (12),  the  trajec¬ 
tories  in  spatial  coordinates  may  be  represented  as  z  =  z(Z(t),  t) .  These  have 
slope 

dz  3_z  3z  dZ  nQ. 

dt  "  3t  3Z  dt  (  ’ 

The  first  term  on  the  right-hand  side  is  given  by  Equation  (15),  anu  the 
second  term  is  the  product  of  (1  +  e)  and  c_  by  Equations  (12)3  and  (14). 
Therefore,  in  spatial  coordinates  the  trajectories  of  constant  strain  have 

speed  as  follows. 

~  =  W  -  f  c  de  +  (1  +  e)  c  .  (20) 

dt  o  JQ  p  P 

Equation  (20)  shows  that  if  c  depends  only  on  e,  these  trajectories  are 
straight  lines  in  spatial  coordinates  as  well  as  material  coordinates.  Several 
lines  of  constant  strain  are  shown  in  Fig.  9.  Note  that  as  penetration 
occurs,  the  higher  strain  levels  are  actually  carried  into  the  target  cavity. 

The  spatial  trajectory  of  any  rod  station  may  be  computed  by  integrating 
Equation  (15)  with  respect  to  time,  Z  being  held  fixed.  All  that  is  needed  is 
the  strain  history  at  a  fixed  rod  station,  and  those  histories  are  available 
eithev  directly  from  the  strain  gage  records  or  by  interpolation  from  Fig.  6. 
Some  trajectories  of  this  kind  are  also  shown  in  Fig.  9.  Note  that  the 
retardation  due  to  axial  compression  is  clearly  visible.  The  trajectory 
labeled  96  mm  corresponds  to  the  initial  length  of  the  recovered  penetrator 
as  determined  by  weighing  in  one  particular  test- 
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Figure  9.  Composite  penetration  diagram  in  laboratory  coordinates 


Finally,  the  spatial  trajectory  of  penetration  is  also  shown  in  Fig.  9. 
This  has  been  drawn  simply  by  translating  the  points  of  the  material  trajectory 
by  the  amount  of  the  free  surface  motion  at  corresponding  times,  ignoring  any 
compression  of  target  material  between  the  penetrator  and  the  free  surface. 
Surface  motion  was  measured  by  Netherwood,  using  a  streak  camera,  and  in  some 
of  his  tests  a  small  plug  approximately  5  mm  in  thickness  was  also  recovered.24 

The  trajectories  of  the  96  mm  station,  the  free  surface,  and  penetration, 
together  with  the  known  plug  thickness,  form  a  remarkably  coherent  picture, 
and  show  with  graphic  clarity  the  need  to  account  for  deformation  in  a  sound 
physical  theory  of  penetration.  That  need  is  evident  in  spite  of  the  uncer¬ 
tainty  in  locating  the  penetration  curve. 

The  highest  strains  measured  all  propagate  with  nearly  the  same  speed. 

In  Fig.  9  that  speed  appears  to  be  slightly  less  than  the  rate  of  penetration 
so  that  it  is  not  clear  from  the  experiments  where  the  erosion  trajectory 
should  be  shown  in  Fig.  8.  A  good  approximation  is  probably  a  straight  line 
drawn  slightly  to  the  left  of  the  highest  strain  shown,  but  with  a  slightly 
greater  slope. 

C .  Experiments  to  Measure  Target  Stress. 

In  a  third  series  of  experiments  Pritchard25'26,27  has  measured  the  stress 
in  the  target  normal  to  the  path  of  penetration  and  directly  on  the  path.  A 
sketch  of  the  experiment  is  shown  in  Fig.  10.  A  manganin  foil  stress  gage  was 

placed  in  the  center  of  a  split  rod,  the  rod  was  inserted  into  a  hole  drilled 

in  the  target,  and  the  alignment  was  adjusted  so  that  the  gage  lay  directly  on 

the  line  of  penetration.  This  experiment  was  also  performed  in  a  4"  light  gas 
gun.  Since  the  manganin  gage  is  subjected  to  a  diverging  strain  field,  it  is 
necessary  to  compensate  by  measuring  strain  independently  and  correcting  the 
pressure  determination.  This  is  done  with  a  Constantin  gage,  which  is  actually 
interleaved  with  the  manganin  gage.  Details  have  been  given  by  Pritchard25'26 

Typical  data  are  shown  in  Fig.  10  for  a  tungsten  alloy  rod  and  a  steel 
target2 '  Stress  histories  begin  with  shock  arrival  and  increase  until  gage 
failure  some  time  before  actual  arrival  of  the  penetrator.  Notice  that  there 
is  no  sharp  jump  in  stress  at  shock  arrival,  at  least  for  these  cases  in  which 
the  penetrator  had  a  hemispherical  nose,  and  that  the  rate  of  increase  is  more 
rapid  for  the  shallower  gage  location.  For  the  two  records  shown,  gages  failed 


24P.  H.  Netherwood ,  Private  Communication ,  1979. 

o  r 

D.  S.  Pritchard ,  "Measurements  of  Dynanic  Stress  and  Strain  Components  in 
Targets  Struck  by  Penetrators,"  ARBRL-MR-03C95 ,  Ballistic  Research  Laboratory, 
APG,  MD,  1981  (AD  A100724). 

260.  S.  Pritchard,  "Piezoresi stive  Gauge  Measurements  in  Deforming  Environments 
During  Penetration, "  Proa.  First  Symp.  on  Gauges  and  Piezoreeistive  Mater¬ 
ials,  Arcachon,  France,  1981. 

27D.  S.  Pritchard  and  G.  E.  Hauver,  Private  Communication  of  Unpublished  Data, 
1982.  The  raw  data  were  taken  by  Pritchard  and  reduced  and  analyzed  by 
Hauver. 
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Figure  10.  Measured  target  stress  vs.  time  at  depths  of  12.7  mm 
and  19.0  mm  below  the  impact  surface.  The  rod  is  a 
tungsten  alloy  and  the  target  is  armor  steel. 
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at  less  than  half  the  time  to  penetrator  arrival  so  that  extrapolation  to 
penetrator  arrival  is  hazardous  at  best.  Nevertheless,  it  is  at  least  inter¬ 
esting  to  note  that  reasonable  extrapolations  give  a  peak  stress  that  is 
roughly  in  the  range  of  50-55  kb  for  both  curves.  Data  have  also  been  reported 
for  the  same  combination  of  materials  shown  in  Fig.  6.  The  general  trend  of 
the  data  is  the  same  as  Fig.  10,  but  the  quality  of  the  data  is  not  as  high. 
(See  References  25  and  26.) 


V.  OTHER  CONSIDERATIONS 

The  experiments  described  above  go  a  long  way  towards  indicating  the 
phenomena  that  should  be  included  in  a  reasonable  theory  of  long  rod  penetra¬ 
tion,  but  there  are  other  important  considerations  as  well.  These  may  be 
grouped  into  three  sets  of  phenomena:  hole  growth  in  the  target,  the  failure 
and  discard  of  eroded  penetrator  material,  and  the  effects  of  bending. 

Neither  the  eroding  rod  model  nor  any  of  the  experiments  described  above 
gives  any  consideration  to  the  dynamics  of  hole  growth,  although  some  attention 
has  been  devoted  to  this  problem  in  the  literature,  e.g.,  Hanagud  and  Ross20 
or  Ravid  and  Bodner.1 4  As  the  penetrator  moves  through  the  target,  target 
material  must  flow  to  the  sides  around  the  penetrator.  At  the  lower  rates  of 
penetration,  the  extent  of  lateral  motion  must  be  at  least  great  enough  to 
admit  passage  of  the  eroded  and  discarded  penetrator  material,  if  any.  At 
higher  rates,  inertia  may  drive  both  target  and  eroded  penetrator  material 
somewhat  away  from  the  sides  of  the  remaining  intact  penetrator,  or  in  the 
case  of  a  noneroding  penetrator,  inertia  may  cause  separation  of  the  flow 
between  target  and  penetrator. 

Ductility  and  fracture  characteristics  also  have  a  significant  effect  on 
penetration.  For  example,  a  very  iuctile  material  may  simply  pile  up  on  the 
target  surface  at  low  impact  velocities  or  actually  turn  inside  out  to  fovm  a 
tube  as  it  penetrates  at  higher  velocities.  (See  Reference  29.)  In  any  case, 
excessive  penetrator  ductility  tends  to  increase  the  contact  area  between 
penetrator  and  target,  thus  providing  a  larger  retarding  force  to  the  penetra¬ 
tor.  In  attempting  to  simulate  numerically  seme  of  Hauver's  data  on  instru¬ 
mented  rods,  Misey30  found  that  it  was  essential  to  include  at  least  a  crude 
fracture  criterion  in  his  two-dimensional  code.  Without  a  fracture  criterion 
the  computed  strains  in  the  penetrator  increased  much  more  rapidly  than  the 
experimental  data  indicated.  This  would  seem  to  indicate  that  a  fracture 
criterion  would  also  be  useful  in  an  engineering  model  of  penetration. 


20  S.  Hanagud  and  B.  Rob8,  " Large  Defamation ,  Deep  Penetration  Theory  for  a 
Canpre88ible  Strain-Hardening  Target  Material , "  AIAA  Journal,  9,  DOS-911, 
1971. 

29 K.  Frank,  Private  Cormunioation  on  teet  result e. 

30e7.  <7.  Misey ,  A.  D.  Gupta ,  and  J.  D.  Wortman ,  " Comparison  of  Penetration  Codes 
for  Strain  Meaourements  in  Kinetic  Energy  Penetrator s ,"  ARBRL-TR-02231 , 
Ballistic  Research  Laboratory,  APG,  MD,  1980  (AD  A0B4391). 
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Finally,  in  any  complete  engineering  theory  of  penetration  some  consider¬ 
ation  must  be  given  co  the  effects  of  bending  as  induced  by  oblique  impact  or 
P  etrator  yaw  at  impact,  and  should  include  lateral  stability  or  bending 
enhancement  due  to  axial  loading.  A  start  was  made  by  Abrahamson  and 
Goodier,31  who  considered  lateral  stability  of  an  elastic/ linear  work  harden¬ 
ing  rod  in  normal  impact  on  a  rigid  target.  However,  they  did  not  consider 
the  effect  of  moments  or  transverse  forces  at  the  impact  end  as  might  be 
induced  by  yaw  and  obliquity,  and  even  more  important  from  the  point  of  view 
of  penetration  mechanics,  they  did  not  include  the  effect  of  penetrator  ero¬ 
sion,  which  in  effect  provides  a  moving  boundary  for  the  application  of  forces 
and  moments.  No  thorough  experimental  study  on  bending  effects  seems  to  have 
been  undertaken,  although  limited  data  exists,  Hauver.1' 


VI.  COMPUTER  CODES 

The  third  approach  to  penetration  mechanics  uses  large  scale  numerical 
simulation  with  the  full  equations  of  continuum  physics.  The  characteristics 
of  some  typical  codes  and  their  limitstions  have  been  described  by  Zukas32  and 
Jonas  and  Zukas33  so  that  only  a  few  comments  are  required  here.  The  ordnance 
regime  still  remains  today  as  extremely  difficult  for  making  accurate  numeri¬ 
cal  predictions,  especially  for  impact  speeds  near  the  limit  velocity.  This 
difficulty  is  commonly  attributed  to  lack  of  adequate  knowledge  concerning  both 
dynamic  material  properties  and  fracture  descriptions  in  the  codes.  (See  Refer¬ 
ence  33  for  example.)  The  determination  of  deformation  and  fracture  properties 
at  high  pressures  and  high  rates  of  strain  requires  exceptional  experimental 
skill  and  ingenuity,  especially  since  the  dynamic  situation  makes  it  all  but 
impossible  to  determine  constitutive  properties  without  imposing  a  priori  con¬ 
stitutive  models  or  assumptions.  There  is  always  some  circularity  in  the 
common  procedure  of  "backing  out"  constitutive  properties  by  choosing  free 
parameters  in  a  given  model  so  as  to  obtain  rough  agreement  between  a  calcula¬ 
tion  and  experimental  data.  On  the  other  hand,  there  often  seems  to  be  little 
other  choice.  All  that  can  be  asked  is  that  critical  judgment  and  restraint 
be  used  in  interpreting  and  accepting  results.  Reference  1  gives  a  thorough 
review  of  current  practice  in  constitutive  modeling  for  codes  and  recommends 
avenues  for  future  improvements. 

As  constitutive  models  become  more  comprehensive,  ever  greater  burdens 
will  be  placed  on  computing  capacity,  which  must  be  offset  either  by  increases 
in  machine  size  and  speed  or  by  increases  in  the  efficiency  of  computing 
algorithms.  Computing  power  should  be  reserved  for  those  parts  of  the  problem 


31 G.  R.  Abrahamson  and  .7.  N.  Goodier ,  "Dynamic  Flexural  Buckling  of  Rods 
within  an  Axial  Compression  Wave,"  J.  Appl.  Mech. ,  33,  241-247 ,  1966. 

32 J.  A.  Zukas,  " Numerical  Simulation  of  Impact  Phenomena"  (Chap  10),  "Three- 
D'imensional  Computer  Codes  for  High  Velocity  Impact  Simulation"  (Chap  11), 
in  Impact  Dynamics,  J.  A.  Zukas,  T.  Nicholas,  H.  F.  Swift,  L.  B.  Greszczuk, 
D.  R.  Curran,  New  York,  Wiley- Inter science,  1982. 

33 G.  H.  Jonas  and  J.  A.  Zukas,  "Mechanics  of  Penetration:  Analysis  and  Experi¬ 
ment,"  Int.  J.  Eng.  Sci.,  16,  879-903,  1978. 
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that  most  need  it,  and  not  wasted  on  regimes  where  liitle  is  happening. 
Fortunately,  rapid  advances  are  being  made  in  both  hardware  and  software  so 
that  advances  in  constitutive  modeling  can  be  expected  to  be  implemented  in 
penetration  calculations. 


VII.  CONCLUSIONS 

The  three  current  approaches  to  penetration  mechanics  -  data  correlation, 
engineering  models,  and  large  scale  computer  codes  -  will  all  continue  to 
evolve  with  increasing  emphasis  being  placed  on  the  role  of  material  properties 
and  the  mechanics  of  material  deformation  and  failure  in  each  case.  In  many 
respects  the  approach  using  engineering  models  is  the  least  developed  of  the 
three,  but  it  holds  the  promise  of  encompassing  virtually  all  of  the  relevant 
physics,  yet  at  far  less  computational  cost  than  the  large-scale  codes.  The 
benefit  comes  from  one -dimensional  rather  than  three-dimensional  numerical 
simulation.  Data  correlations  will  still  be  indispensable  for  organizing 
large  volumes  of  ballistic  data,  but  engineering  models  can  be  expected  to 
suggest  improved  functional  forms  for  non dimensional  analysis.  The  engineer¬ 
ing  models  themselves  should  play  their  biggest  role  in  aiding  designers  and 
systems  analysts.  The  large  codes  will  remain  as  the  indispensable  tool  for 
examining  tha  details  of  interaction,  and  as  constitutive  models  improve,  may 
well  permit  actual  material  design  for  ballistic  applications. 
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